When do tracer particles dominate the Lyapunov spectrum? 
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Dynamical instability is studied in a deterministic dynamical system of Hamiltonian type com- 
posed of a tracer particle in a fluid of many particles. The tracer and fluid particles are hard balls 
(disks, in two dimensions, or spheres, in three dimensions) undergoing elastic collisions. The dy- 
namical instability is characterized by the spectrum of Lyapunov exponents. The tracer particle is 
shown to dominate the Lyapunov spectrum in the neighborhoods of two limiting cases: the Lorentz- 
gas limit in which the tracer particle is much lighter than the fluid particles and the Rayleigh-flight 
limit in which the fluid particles have a vanishing radius and form an ideal gas. In both limits, a gap 
appears in the Lyapunov spectrum between the few largest Lyapunov exponents associated with the 
tracer and the rest of the Lyapunov spectrum. 



During the last decade, dynamical instability and chaos in systems of interacting particles has become a problem 
of major preoccupation in statistical mechanics. Many systems have been shown to present sensitivity to initial 
conditions characterized by positive Lyapunov exponents [§],||. Methods from kinetic theory have been developed 



transport properties have also been investigated |1C|-|13|| . 

In mixtures of a one-component fluid with a very low concentration of identical tracer particles the largest Lyapunov 
exponent may either be virtually identical to the largest Lyapunov exponent of the pure fluid (with just a slight 
perturbation caused by the tracer particles), or it may be larger, due to the dynamical properties of the tracer 
particles. In the latter case one may say that the largest Lyapunov exponent is dominated by the tracer particles. 

The purpose of the present paper is to describe the regimes in which the tracer particles dominate the dynamical 
instability of the fluid system. If the mass of the tracer particle is not much larger than that of the fluid particles, and 
in addition the mean free path of both bath and tracer particles are at least comparable in length to their respective 
radii, the characteristic Lyapunov exponents for fluid and tracer dynamics respectively are roughly of the order of 
their respective collision frequencies. So in this case one may conclude that the tracer particle will dominate if its 
collision frequency is sufficiently larger than that of the fluid particles. 

One obvious way to reach this goal is by making the mass ratio M/m of tracer mass over fluid particle mass very 
small. The limit where this ratio goes to zero corresponds to the Lorentz gas, for which the Lyapunov exponents of 
the tracer particle were derived several years ago f|||H ■ 

Another common situation is where the radius of the tracer particles is much larger than that of the fluid particles, 
e.g., when one satisfies the conditions for Brownian motion, where also the mass of the tracer particle is much larger 
than that of the fluid particles. In this case the collision frequency of the tracer particles is much higher than that of 



I. INTRODUCTION 




11 



,?]. The relationships to 



f 



the fluid particles, but at the same time the velocity changes of the tracer particle in collisions with a bath particle 
are much smaller than those of the bath particles themselves. In addition the curvature of the Brownian particle, 
due to its large radius, is much smaller than that of a fluid particle and therefore the diverging effect (crucial for a 
positive Lyapunov exponent) of a Brownian-fluid collision is much smaller than that of a fluid-fluid collision. As a 
result of this the largest Lyapunov exponent usually is determined by the fluid particles, as noted both by Louis and 
Gaspard jlij and by Nasser and Dorfman jt5). However, there is one noteworthy exception to this, corresponding to 
the so-called Rayleigh-flight limit, where at fixed Brownian mass and radius the radius of the fluid particles is sent 
to zero. In the extreme limit of a Brownian gas surrounded by an ideal gas it is obvious that the Brownian particle 
has to dominate, as the Lyapunov exponents of the ideal gas are strictly zero. The goal of the present paper is thus 
to study these cases of dominance of the Lyapunov spectrum by the tracer particle. 

The plan of the paper is the following. In Se c. O , the problem of dynamical instability in a system of hard balls 
of different masses and radii is posed. In Sec. |ni| we consider the Lorentz-gas limit. In Sec. [Fv| we consider the 
Rayleigh-flight limit. Conclusions are drawn in Sec. fv[ 



In this paper we will consider a system composed of many fluid particles and one tracer particle. From a general 
viewpoint, this system can be considered as a binary mixture of Nf fluid particles, which we take to be hard spheres 
of radius a and mass m, with JVt tracer particles, which will likewise be hard spheres, with radius A and mass M. All 
the N = N{ + N t hard balls move in a rectangular domain of finite extension £1 with periodic boundary conditions. 

The motion of the hard balls is composed of free flights between binary collisions which are elastic and instantaneous. 
Energy and the total linear momenta are conserved. 

For a system of hard balls of radii {ai}fL 1 and of masses {rrii}fL 1 , the equations of motion are given as follows in 

terms of the positions and velocities {rj \t n ),v^ (i«.)}£Li an d { r i v i(in)}£Li> respectively before and after 

the collision at time t n : 



II. THE DYNAMICAL SYSTEM AND ITS INSTABILITY 



A. The dynamics 



1. Free flight between binary collisions: 




(1) 



2. 



Binary collision: 



n 





(2) 
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with the unit vector joining the centers of the i and j balls at the collision given by 

_(±) __(±) 

i j 

£ij — ; 

cii + aj 



(4) 



and the relative velocity vector 



v|7 ) = v<->-v<"> (5) 



B. The linearized dynamics 

The dynamical instability of this system is characterized by the rates of exponential growth of infinitesimal per- 
turbations on the positions and velocities of the particles: <5X = {ftr^tfvi}^. The rates are called the Lyapunov 
exponents 

Depending on the initial perturbation <5X(0) we may have as many different Lyapunov exponents as there are inde- 
pendent directions in phase space. Since the dynamics of the present hard-ball system is of Hamiltonian (symplectic) 
character the Lyapunov exponents obey a pairing rule: If Xi is a Lyapunov exponent, then — is also a Lyapunov 
exponent. The set of exponents form the so-called Lyapunov spectrum. The Lyapunov exponents associated with the 
directions perpendicular to the energy and momenta shells vanish. In the present system, which has no fixed points 
except at zero energy, this also holds for the pair mates associated with the directions of time and center-of-mass 
translation. Accordingly, 2 + 2d Lyapunov exponents vanish. 

The problem of the dynamical instability of a hard-ball system was formulated in the seventies by Sinai ||l6|| who 
was inspired by the pioneering work of Krylov in the forties ]T^ ]. Independently, Erpenbeck and Wood carried out 
numerical investigations, also in the seventies []l8|| . The systematic calculation of the Lyapunov spectrum in hard-ball 
systems has been developed in the nineties. Using the method of Gaspard and Dorfman jl9| one can derive the 
following linearized equations in terms of the infinitesimal perturbations before and after each collision | fl4[ : 

1. Free flight between binary collisions: 



2. Binary collision: 
gj+) - xj-) _ o m i (r Sr^) c 

Sr^ = <5rH + 2 ^Me l? • Sr^h e ■ (8) 
<5r[, + ^ = i5rj^ for k^i,j 
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6vtf = 5vH - «vH 



(9) 



(10) 



(11) 
(12) 



C. The kinetic properties of the thermodynamic equilibrium state 



We require that the center of mass is at rest and we are interested in the properties of the equilibrium thermodynamic 
state at fixed temperature T. Accordingly, the total linear momenta vanish: Ptot — J2iLi miVi — 0, while the total 
energy is given by 



N 1 rl 
E = "£-mrf = ~(N-l)k B T 



(13) 



where fee is Boltzmann's constant. For hard-ball systems, the motions at different temperatures are equivalent up to 
a rescaling of time. In the sequel, the temperature is thus fixed at the value T = fcg 1 . 

In the fluid phase at low enough density, the system is supposed to be ergodic on each energy-momenta shell, which 
defines the equilibrium states. The mean velocity^ of each particle is defined as 



(IN 



(14) 



At equilibrium, the mean velocities are determined by the temperature and the mass of the particles. For large particle 
number, they are given, to leading order in 1/JV, by 



d = 2: 
d = 3 : 



2m,i 
8k^T 



(15) 
(16) 



where d is the space dimension. 



1 The mean velocity should be well distinguished from the root-mean-square velocity, defined as -J (vf). 



4 



The mean relative velocities between the particles entering a binary collision will also be of importance in the sequel. 
They are defined by 

Vij = (\\ Vi - Vj \\) (17) 

and are given by Eqs. (|l5|)-(|l6|) with the mass replaced by the relative mass 

_ miirij 
^ rrii + rrij 

If the system is sufficiently dilute the collision frequencies of the fluid and tracer particles can be evaluated by 
supposing that each particle has a cross-section for collision with each one of the other types of particles: 

d = 2 : (J i:j =2 (a, + aj) (19) 
d = 3 : Oij = 7r (a» + aj) 2 (20) 

If as, a{ t = ati, and a t t denote the cross-sections for fluid-fluid, fluid-tracer, tracer-fluid, and tracer-tracer collisions, 
the collision frequencies of the tracer and fluid particles to leading order in the fluid density are given by 

N f N t - 1 

vt = m + vtt — -jj- c t f v t f H ^ — a tt v u (21) 

Vf = m + "ft - ere vg + — att «ft (22) 

where vtf — Uft and the extension parameter is the area respectively the volume of the system: 

d = 2 : fi = L x L y (23) 
d = 3 : Q = L X L y L z (24) 

We notice that each tracer particle may collide on the N t — 1 other tracer particles, which explains the presence of 
— 1 in the tracer-tracer collision frequency v t t- The term —1 is important in systems with a low number N t of tracer 
particles. In contrast, it can be neglected if the number of particles is large as it is the case for the fluid particles. 

D. Simulations and system preparation 

In order to test our theoretical results we performed several MD simulations in which we computed the Lyapunov 
spectra of hard ball systems containing a tracer component. In all these simulations, we consider a fluid with a single 
tracer particle, so Nt = 1. 

To initialize a simulation we locate the fluid particles on the lattice points of a crystal lattice. In d = 2, the initial 
positions form a triangular lattice with M x M y rectangular cells of two disks each. In d = 3, the initial positions form 
a face-centered cubic (FCC) lattice with M x M y M z cubic cells of four spheres each. The sizes of the cells are fixed in 
order for the fluid to have a fixed density n in absence of the tracer particle. In d = 2, the domain is rectangular of 
sizes 

d = 2: L X = M X (J^j 2 and L y = M y f^pj (25) 
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so that its area is O = L x L y = 2M x M y /n. In d = 3, the domain is also rectangular but of sizes 

d = 3 : L X = M X (J^J , L y = M y (J^J * , and L z = M z (J^J ' (26) 

so that its volume is O — L x L y L z — 4M x M y M z /n. 

Thereafter, the tracer particle is placed in the middle of the crystal configuration under removal of all the fluid 
particles that would overlap with the tracer particle. The number of fluid particles is thus given approximately by 

N { ~ n (fi - fi t ) (27) 

where 

d = 2 : Q t = 7T A 2 (28) 

4-7T 

d = 3: n t = — A 3 (29) 

As a result of this procedure the density of the fluid particles in the free volume is still given to an excellent ap- 
proximation by n, even if the tracer particle is so large that it occupies an appreciable fraction of the total volume 
0. 

Accordingly, the collision frequency of the tracer particle is 

v t =n otf v ti (30) 

and its mean free path is 

4 = ^-!*— (31) 
v t n CT t f w t f 

Table I gives these quantities in d — 2 and = 3. 



Table I. Characteristic quantities for the motion of a tracer in a dilute fluid. 



dimension 
mean velocity 

collision frequency 
mean free path 



d=2 



2M 



p. = a*. 1 r 

Ct v t — 2(a+A)n\j r 



m+M 



d = 3 



Sk B T 
M 



_ [&k 

~ V ~ 

/ . ,,x2 / 8k B T(m+M) 

t — ~ i r 

~~ vt — ir(a+i) a n y r 



m+M 



Finally, the initial velocities of the particles are drawn with a random number generator from a Maxwellian distri- 
bution with the appropriate mass. In practice we choose fcflT 1 = 1 as well as m = 1. 



G 



E. Discussion of the dominance of the dynamical instability by the tracer 



Since there is a single tracer particle it only collides with the fluid particles. However, each fluid particle may 
collide both with other fluid particles and with the tracer particle. Accordingly, there are two types of collisions: the 
fluid-fluid collisions and the tracer-fluid collisions. 

Since both types of collisions happen between particles with convex surfaces, we expect that both of them will 
contribute to the dynamical instability. As discussed in the introduction, the strongest instability, characterized by 
the largest Lyapunov exponent, could be dominated either by the fluid- fluid collisions, or by the tracer-fluid collisions. 
To decide which of these possibilities is realized for a given choice of parameters we can proceed in the following 
way: We calculate both the fluid-fluid and the tracer-fluid maximal Lyapunov exponent under the assumption that 
indeed fluid- fluid collisions repectively tracer- fluid collisions are dominant and then compare the results. In the great 
majority of cases indeed the larger of the two calculated exponents gives an excellent approximation to the actual 
value of the maximal Lyapunov exponent. 

The first calculation concerns the fluid-fluid collisions. In the absence of tracer, such collisions give the maximal 
Lyapunov exponent 



Af ~ cu(N{) Vf In 



a(N{) 



with V{ ~ n as vs (32) 



4na d 

where, for disks in d — 2, Van Zon et al. |6j have shown that the prefactor u>{N) is well fitted by the expression 

u w - 4 - 311 - ( 33 ) 

although the actual asymptotic behavior for large N is as I/(logA) 2 jjj. Further, a(N{) is of order unity, depends 
only weakly on Nf and approaches a constant for Nf — > oo. 

On the other hand, we may consider the motion of the tracer particle undergoing kicks by independent fluid particles. 
The equations of motion are given by Eqs. (|)-(|) in which r i=1 = R, v i= i = V. We denote by (rJ, \v£ } ) 
and (R„ ,V„ the position and velocity of the tracer particle before repectively after the n th collision, which 
occurs at the time t n . Since the trajectory is continuous in position, R4 = R« = R„. On the other hand, 
v£ + ' ) = V^ + \ = V„ . We have the following iteration for the motion itself: 

R n = R n _ 1 +r n V n _i (34) 
2 m 

V n = V„_i — — [e„ • (V„_i - v„_i)] e„ (35) 

m + M 

tn=tn-l+T n (36) 

where r n is the time interval between the (n — l) th and the n th collision, while 

R„ 



A 



(37) 



is the unit impact vector at the n th collision. 

An infinitesimal perturbation on the motion is ruled by Eqs. (|7])-([l2]) in which 5rj = i = 5R, 5v,-i = 5V and 

$Tj = 5Vj = for j =/= 1 because, if indeed the tracer particle dominates the maximal Lyapunov exponent, the 
perturbations of the fluid particle positions and velocities will be negligibly small compared to those of the tracer 
particle. For this argument to hold it is important that recollisions of a fluid particle with the tracer particle are either 
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rare, or, if they are not, are still dominated by the per turbations of the tracer particle. Especially when M 3> m this 
is a subtle point, because, as can be seen from (fj])-(]l^), the perturbations of a fluid particle right after a collision are 
comparable to those of the tracer particle. So if a recollision is not unlikely, it has to occur typically after a time that 
is long, compared to the Lyapunov time of the tracer particle. The iteration for an infinitesimal perturbation on the 
motion is then given by 



free flight: 






+ r n (TV n _ 


binary collision: 




= «h4 _) - 


2 m 
m + M 








2m 
m + M 


with SV n = M +) = <5vi+i 


and 







(e n ■ SR^ e„ 

{(e„ • <5V„_i) e„ + [<5e„ • (V n _i - v n _i)] e„ 
+ [e„ • (V„_i - v n _i)] 5e n } 



A + a 



6R ( -*> - (V»_i - v„_i) 



(-) 



(V n _i - v n _i) 



(38) 
(39) 

(40) 
(41) 



In the equations above, the perturbation SR n ' corresponds to the first impact time of the two neighboring trajectories 
and the perturbation SR n +S> to the later impact time. They are related as 



<5R« = 0H4-) + (V n _i - V„) Sr n 



with the time lag at collision 



0~T n = - 



e n ■ SR. 



C-) 



£n ' (V n -1 ~ V n -l) 



(42) 



(43) 



These equations can be analyzed thanks to a few observations. First of all, the time r„ between two successive 
collisions has an average value given as the inverse of the collision frequency of the tracer particle: 



(rn) = ~ 

Furthermore, we find the mean relative velocity as 

(l|V„_i - v n _i||) = vti 



(44) 



(45) 



We may expect a dominance of the resulting Lyapunov exponents for the tracer particle over those resulting from 
the fluid-fluid collisions in two cases to be discussed below: in the Lorentz-gas limit M — > with M <C m, and in 
the Ray leigh- flight limit a — > with a -C A and na d -C 1. The calculation of these Lyapunov exponents will be the 
subject of the next two sections. 



III. THE LORENTZ-GAS LIMIT 



In the limit M — > 0, with M <C m, the tracer particle is much faster than the fluid particles. Accordingly, the tracer 
moves through a fluid which is essentially at rest. This system is referred to as a Lorentz gas pO|,pT|. The collision 
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frequency of the tracer particle will thus be larger than the collision frequency of the fluid particles. Therefore, the 
perturbations on the coordinates of the tracer particle will grow faster than the perturbations on the fluid particles. 
Actually the value of the maximal Lyapunov exponent of such a fluid can be predicted in this case thanks to the 
work by Van Beijeren, Dorfman, and Latz H,^). 



A. The two-dimensional case 



In d = 2, the Lorentz gas has a single positive Lyapunov exponent as a consequence of chaos and energy conservation. 
Therefore, we expect that the maximal Lyapunov exponent is the positive Lyapunov exponent of the Lorentz gas when 
M -C 771, while the next Lyapunov exponents remain at the level of the fluid Lyapunov exponent (|3^). For a dilute 
system with M <C m, the maximal Lyapunov exponent is therefore given by 



d = 2 



Ax ~ 2 (a + A) n 



irk B T(m + M) 



2mM 



In 



,i-c 



2n(a + A) 2 



(46) 



with Euler's constant C = 0.5772156649... g|. 

This behavior is well confirmed by numerical computation. Figure [l] compares Lyapunov spectra for M 3> m and 
for M <C m at a fixed ratio of a and A. We observe that the maximal Lyapunov exponent is separated from the rest 
of the spectrum for M <C m. This implies that with increasing m/M a gap opens up in the Lyapunov spectrum. 

Figure || shows that in the limit M — ► the collision frequency of the tracer particle increases as predicted by Eq. 
@ in the limit M -> 0. 

Figure |^ depicts the five largest Lyapunov exponents as a function of the mass of the tracer particle showing that, 
indeed, the Lyapunov spectrum is dominated by the tracer particle as soon as M <C m. The maximal Lyapunov 
exponent undergoes a cross-over from the fluid value (^) for M > m to the Lorentz-gas value (|46) for M <C m. In 
the regime M <C m, the second Lyapunov exponent tends to slightly increase up to the fluid value (32|) A2 — Af < Ai. 
A similar behavior is seen for the next Lyapunov exponents. 

Figure ^ shows the dependence of the spectrum on the radius A of the tracer particle. As the tracer becomes larger 
and larger in a fluid of fixed density, the room for its motion between the fluid particles becomes relatively smaller 
and smaller. A cage effect occurs for a large tracer as though the effective density of the Lorentz gas would increase 
(indeed the motion of a mobile particle of radius A among fixed point scatterers is completely equivalent to that of a 
mobile point particle among fixed scatterers of radius A that are allowed to overlap each other. In the literature this 
has been treated as the Lorentz gas with overlapping scatterers). In this case, the assumption it [A 3> 1, on which the 
demonstration of (^6|) is based, breaks down, which explains the discrepancy observed in Fig. Q for 4»a between 
the numerical results and the prediction of Eq. (p6j) . 



B. The three-dimensional case 



In d — 3, the Lorentz gas has two positive Lyapunov exponents as a consequence of chaos and energy conservation. 
In this case, we thus expect that the two largest Lyapunov exponents are the positive Lyapunov exponents of the 
Lorentz gas when M <C m, while the next Lyapunov exponents remain at the levels of the fluid Lyapunov exponents 
(|32"|). For a dilute system with M <C m, the two largest Lyapunov exponents are therefore given by || 



d= 3 : 

Ai ~ 7T (a + Af n 



8fc B T(m + M) 



In 



4e" 



nn[a 



A) 3 



(47) 
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A 2 ~ 7T (a + A) 2 n 



8k B T(m + M) 
irmM 



In 



r-C 



mr(a + A) 3 



(48) 



The next Lyapunov exponent should remain at the fluid value (|3^): A3 ~ Af < A2 < Ai. 

Here again, this behavior is well confirmed by numerical computation. Figure |B| depicts a Lyapunov spectrum for 
M<Cm where we observe that, indeed, two Lyapunov exponents are separated from the rest of the spectrum, again 
with the formation of a gap as m/M increases. Figure [| depicts the five largest Lyapunov exponents as a function 
of the mass of the tracer particle showing that the Lyapunov spectrum is dominated by the tracer particle as soon 
as M < to. The two largest Lyapunov exponents undergo a cross-over from the fluid values when M > m to the 
Lorentz-gas values (^) and ( ^8[ ) for M <C to. In the regime M -C to, the third and next Lyapunov exponents tend 
to slightly increase up to the values of the first and next Lyapunov exponents of the previous regime M»m. 



IV. THE RAYLEIGH-FLIGHT LIMIT 



Another regime in which we may expect that the tracer particle dominates the Lyapunov spectrum is the one near 
the limit where the radius of the fluid particles vanishes: a — ► 0. This regime can be characterized by the conditions 
A ^> a and na d <C 1. In the limit, the fluid particles have no collisions between each other, but only with the 
tracer particle. This is referred to as the Rayleigh-flight limit and is known to present diffusive motion 23 1. In 
the limit a — > the fluid Lyapunov exponents ([52]) vanish. The tracer particle undergoes elastic, diverging collisions 
with the fluid. Since the tracer dynamics has d degrees of freedom, and the presence of the fluid particles precludes 
symmetry transformations of the tracer coordinates alone, we expect the existence of d positive Lyapunov exponents. 
In the Rayleigh-flight regime, we may thus expect that the dynamical instability of the fluid is dominated by the 
tracer particle, which is confirmed by the analysis given below. We want to remark here that the Rayleigh-flight limit 
provides an example of a system where the Lyapunov exponents remain well-defined in the infinite system limit. In 
fact taking this limit simplifies the analysis, because it largely limits the possibilities of recollisions. 



A. Dynamical instability in the Rayleigh flight 

Before considering the problem of the full fluid with a ^ 0, let us consider the dynamical instability of the tracer 
particle in an infinite domain filled with an ideal gas of fluid particles without mutual interaction. The fluid particles 
of mass to are coming from infinity with velocities distributed according to a Maxwell-Boltzmann distribution with 
temperature T and have a uniform spatial distribution with density n. In the limit a — > the fluid particles have an 
infinite mean free path. The only possible collisions occur with the tracer particle. We notice that a fluid particle may 
collide more than once with the tracer particle. However, recollisions are rare if either of the conditions M/m 3> 1 or 
lt/(A + a) 3> 1 are satisfied. It will turn out that in order for the tracer particle to be dominant, at least one of these 
conditions has to hold and therefore we will neglect recollisions in the sequel. 

At each collision, the velocity v = v„_i of the fluid particle is a random vector of Maxwell-Boltzmann distribution 

. . / TO \ 2 / TOV 2 \ 

P W = (ster) '(m) (49) 

Since the spatial distribution of the fluid particles is uniform the impact positions are distributed uniformly over the 
cross-section. If <j> = (j> n denotes the angle between the impact unit vector e n with the relative velocity v n _i — V n _i 
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the uniform distribution of the impact positions implies that the collision angle is distributed according to 

d=2: P{(f)) = cos <f) (51) 
d = 3: P((/>) = sin20 (52) 

for (j> G [0, . At low density and for M ^> to, the successive collisions undergone by the tracer particle occur at 
random time intervals r = r„ which, for given speed V — ||V„_i|| of the tracer particle are distributed according to 
an exponential probability distribution of mean intercollision time (t(V)) = l/u t (V): 

P(r) = u t (V) exp[-u t (V)r] (53) 

In the limit M/m — > oo the velocity dependence of the collision frequency disappears. 

If the recollisions are neglected, the parameters (v„_i,r„,e„) of the successive collisions are independent random 
variables, though for finite mass ratio their distribution does depend on the tracer speed V. Therefore, the tracer 
particle follows a random process which can be numerically simulated by a Monte-Carlo method. The growth factors of 
infinitesimal perturbations on the position and velocity of the tracer particle also follow a stationary random process. 

Two regimes can be distinguished, depending on the value of the parameter 

to 1 j to it 

1 ~ M + m n{A + a) d ~ V M + to A + a ^ ' 

This parameter is always smaller than the ratio of the mean free path i t of the tracer particle to the collisional radius 
A + a since to/ (to + M) < 1. 

For 7> 1 the dynamics in tangent space is very similar to that of a Lorentz gas. Typically the position perturbation 
of the tracer particle when entering a collision with a fluid particle, can be approximated by the product of its velocity 
perturbation and the free flight time since the previous collision. Since it 3> A + a in this case, the dynamical 
instability is essentially dominated by the large distance to the next collision that amplifies the perturbation on initial 
conditions by the ratio it /(A + a) (see Fig. Ma). Therefore, the maximal Lyapunov exponent can be calculated again 
by just considering the dynamics from one collision through the next one and averaging over the parameters of such 
an event. In fact, the Lorentz gas regime may be considered as a special subset of the class of all systems meeting the 
requirement 7> 1. This condition alone, together with the requirement that the resulting maximal tracer Lyapunov 
exponent exceeds the maximal fluid Lyapunov exponent, is sufficient to have dominance of the tracer particle. 

For 7 « 1 typically also the relative changes of the perturbations in a collision are much smaller than unity. A 
possible typical situation is illustrated in Fig. [7]b in this case. Therefore, we may describe the dynamics in tangent 
space to a good approximation by a Fokker-Planck equation. Even though we have not found an analytical solution 
for this, we can infer several important properties from it, especially we can find out how the maximal Lyapunov 
exponent scales with the parameter 7. 



(a) The regime 7>1 



In this regime, we can use Krylov's argument |l7| ] to estimate the maximal positive Lyapunov exponent as follows. 
If 7 S> 1, Eq. ( |4l| ) shows that a perturbation on the impact unit vector can be estimated as 

Se n ~ S9 n ^i (55) 
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where S9 n -i is a perturbation on the angle of the tracer velocity. If we substitute this estimation into Eq. 
can conclude that the perturbation on the tracer velocity angle should grow on average as 



Wtf 



AI + m A + a v t 



we 



(56) 



in a collision. Accordingly, the maximal Lyapunov exponent behaves as A t ~ vt In 7. 

The Monte-Carlo simulation confirms this expectation. Using the last expression of Table I in Eq. 
that the numerical results for the maximal Lyapunov exponent of the tracer are well represented by 



7 > 1 



A t ~ v t In 



at 



M + m n(A + a) d 



6|), we find 



(57) 



with the constant at = exp(l — |C)/(27r 1 / 2 ) = 0.322 . . ., in the limit M ^> m for d = 2. This agrees well with our 
simulation results, which give a fitting value for at of 0.33 ± 0.05. In the Lorentz-gas limit M <C m, the quantity at 
also tends to a constant value, given by a t = e 1 " c /2 = 0.763 ... for d = 2, as may be seen from Eq. (Eft) . The second 



positive Lyapunov exponent X' t turns out to have a different dependence on A 
theoretically yet. 



a which we have not investigated 



(b) The regime 7 <C 1 



In this regime, the dynamical instability is weaker. The perturbation on average only grows by a factor slightly 
larger than unity at each collision. From the Lyapunov instability, together with the relationship -^.SR — SV, it is clear 
that, in the mean, the velocity and position perturbations of the Brownian particle are related through <5V ~ A t <5R. 
Due to the Brownian fluctuations in the system, SV will in fact constantly fluctuate around the value A t 5R, both 
in magnitude and in direction. To characterize these fluctuations we develop a continuous-time description of the 
dynamical instability similar to the one developed by Van Zon in another context pif] . We introduce the variables x 
and y, defined through 



5V-x6R 6V ± 

(59) 



||<JR|| pR|| 

The variables x and y = ||y|| describe the ratio's of velocity to position perturbations of the tracer particle, for the 
components parallel and perpendicular to <5R respectively. 
They allow us to obtain the maximal Lyapunov exponent as 

A t = limim''^" 



T ||<5R(0)|| 



1 



T 



= T-Soo T dt X{t) (60) 







because 
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x=j t ln\\SIL\\ (61) 

and, again, ^<5R = SV. Accordingly, the maximal Lyapunov exponent is given by the time average of the variable 
x. By ergodicity, this time average is equal to the statistical average over the stationary probability distribution of x: 

At = (x) (62) 

During each free flight between two collisions, the variables x and y have the form 

ar + b 



q T 2 + 2bT + C ^ 

±Vac — b 2 ,„ A . 
y = ar 2 + 26r + c (64) 



with r — t — t n -\ and 



a=(6V n - 1 f, b = SB™ 1 -5V n - 1 , and c=(SR ( +\) 2 (65) 



The ratio x/y is a linear function of time so that 

K=± (66) 
dt \yj y K ' 

is a constant for the time evolution of x and y during the free flights. Without collisions, the variables x and y thus 
satisfy equations of motion 

ifroo = -x 2 + y 2 = +y 2 d y K 

y flcc = -2xy = ~y 2 d x K (67) 

However, the variables x and y undergo a jump at each collision of the tracer particle with the fluid particles. 
Accordingly, the time evolution of these quantities is ruled by the following coupled stochastic equations 

x = -x 2 + y 2 + i coll (68) 
y = -2xy + y co n (69) 



with 



■ sv coU 

X coll = ^2 ( 70 ) 

5V co ll - p (p ■ JVcoll) . . 

yc ° ii= m\ ' ( } 



where wc introduced the unit vector 



5R 

ps W (72) 



which is orthogonal to the unit vector cr = y/y. 
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The jumps in the perturbation on the velocity are determined by Eq. (|40|) in the following form 

+00 

<5V coU = J2 {SVn-SVn-l)S(t-t n ) 

TL— — CO 

Accordingly, the collisional contributions (|70|)-(|7l|) take the forms 

icon = £ n 5(t - t n ) 



(73) 



ycoii 



+00 

E 



In *(*-*n) 



(74) 



(75) 



To simplify our analysis we assume in the sequel that the mass ratio m/(M + m) is very smalljj As a result the 
relative velocity between tracer and fluid particles may be approximated by minus the velocity of the fluid particle, 
the collision frequency of the tracer particle becomes velocity independent and also the stationary distribution of the 
variables x and y becomes independent of the velocity of the tracer particle. The expressions for the jumps £ n and 
rj n simplify to leading order in M ™ m <C 1. In particular, in Eq. ( ff"o| ) the terms proportional to 5V may be neglected. 
In addition, in relative velocities the contributions of the tracer particle may be ignored. Using Eq. (flO|), we obtain, 
with these simplifications: 



2m 



(to + M)(a + A) 
2m 

(to + M)(a + A) 



( e n ' Pn 



On • Pnf 



Pn 



e rt ' Pn 



e « + (Pn ' U ") e n ~ (En ' P„) U rl 



rjn T % 



with the unit vector 



Vn-l 
|v n -l| 



(76) 
(77) 

(78) 



and with t„ a unit vector orthogonal to p n . Since the successive jumps may be assumed to be statistically independent 
random variables, the time evolution in tangent space of an ensemble of tracer particles may be described by a 
Boltzmann equation for the distribution of x and y, of the form 



/+00 />oo rl 

dU dr, drP(Z,ri,T)lf(x-£,y-nT,t)-f(x,y,t)] 
-00 Jo J 



(79) 



Here "P(£, 77, t) is the distribution function for jumps £ and tjt in x respectively y as result of a collision of the tracer 
particle with a fluid particle. The form of this distribution may be obtained from Eqs. (p0|)-(|4l|), combined with 
the Maxwell distribution for the fluid particle velocity v„_i and the distribution of the impact vector de n for given 



2 Even if this is not the case our main results, notably the way in which the maximal Lyapunov scales with 7, remain valid. 
But the specific analysis is more complicated, because the simplifications made at this point do not apply. 
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V„_i — v„_i. We used here that the terms in Eq. ( [ICf ) proportional to (e„ • (5V„_i) e„ may be neglected, as result of 
which V(£, 1], r) becomes independent of x and y. This is justified because of our condition 7 <C 1, as can be seen on 
hindsight by comparing this term to the remaining terms in jio]), approximating (5R by <5V/At in the latter [see Eqs. 

^he P r i me on the integral over r indicates an integration over the d— 1-dimensional unit sphere orthogonal 
to p. One readily sees that in fact "P(£, rj, t) as function of r is distributed uniformly over this unit sphere. Moreover, 
this distribution has long tails due to so-called grazing collisions, i.e., collisions with small scattering angles. As one 
sees from the second term in ( pl| ) these may give rise to arbitrarily large jumps in x and y of the order of 

£n , Vn ~ (80) 
COS 0„ 

As a consequence of the distribution (pl|), in d = 2, or (p2), in d = 3 the tails of the probability density behave as 



"t 7' 



for |£| 



P(e,^r)~^ JI5 . S (81) 



v t 1 



for |^| — > 00 



Therefore the first moments of this distribution exist. However, the slow decay of the distribution for large jumps 
leads to logarithmic divergences of the second moments of V . 

When the jumps in x and y at collisions are mostly small compared to the typical magnitudes of x and y themselves, 
it seems appropriate approximating the Boltzmann equation by a Fokker-Planck equation. For this to be the case we 
again need the condition 7 < 1. But in addition, the approximation of the Boltzmann equation ([79]) by a Fokker- 
Planck equation requires the existence of the second moments of the distribution 7 5 (^, 77, t). As we just noticed, these 
moments do not exist. But their divergence is only logarithmic and therefore the process, to leading order in 7, can 
still be described by a Fokker-Planck equation. The typical scale for the jumps in x and y is of order £, 77 ~ i/ t 7 with 
7 <C 1. Consequently, if we introduce a cut-off in the calculation of the second moment at values of £ and r\ satisfying 

Ccut-off, ?7cut-off ~ v t (82) 

with some 5 > 0, most of the distribution V falls within the cut-offs. We will justify this cut-off below by an argument 
of self-consistency, and at the same time determine the appropriate value of S. Utilizing the cut-off, we obtain the 
following Fokker-Planck equation: 

dtf + d x [(-x 2 +y 2 + + 8 y [{-2xy)f] = D„ d 2 J + D x± d 2 y f (83) 

where the drift in x is given by the first moment of the distribution V as 

H = * (0 (84) 

and the "diffusion coefficients" are defined as 

% = ^t<(£-<£» 2 > (85) 

D„ = ^ * (V 2 ) (86) 

where d is the space dimension. There is no drift in y because (77) = by symmetry. For the same reason, the 
coefficient of cross-diffusion in x and y vanishes. 

Furthermore, in the limit M/m — > 00 we are considering, one has (V co n) = by isotropy, irrespective of V and R. 
Therefore also (<5V co n) and, according to (|70|) and (|74|), (x C ou) and (£) vanish. A non-zero fj,\\ may be obtained by 
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placing, the Brownian particle in a harmonic well, a case we are investigating presently. ^From the expressions ( |76| ) 
and (|77| ) and the cut-offs (|82|) the diffusion coefficients may be estimated to behave as 



Aul = ~K7~T~ TT *t (»7 2 ) ~ ^ 7 2 In -, 
2(d-l) 7 

Note that equation (|3|) has solutions that are independent of t. It is precisely those solutions that we are interested 
in. 

Continuing with hard spheres we note that, with vanishing first moments, the Fokker-Planck equation may be 
rescaled by applying the scalings i^ikj, y^ycs, and t — > t This allows us to extract the dependence of 

the variables on the expression k = v% ln(l/7). The only remaining parameter that might change on varying the 
ratio's of masses and diameters is the ratio D±±/ D^. But even this will change just very slightly, and in fact only 
does so because the cut-offs on £ and rj in determining the second moments introduce a 7-dependence which works 
out somewhat differently for Dyn and D±±. In the limit k — > 00 this ratio has a well-defined limit. 

From this scaling we may conclude immediately that (x), and hence the maximal Lyapunov exponent, scales as 
k 1 / 3 and we finally obtain that 

A t - (x) ~ u t (V In ^ 3 (89) 

with the parameter 7 defined in (jH^). 

A discussion is now in order to justify the introduction of the cut-off. As mentioned already, the second moments of 
the distribution V(£,ti, t) diverge logarithmically. This implies that in principle we should stick with the Boltzmann 
equation ([79]). Still, the typical scale for the jumps in x and y is of order 7, whereas the typical scale for variations 
of x and y is expected to be at least approximately of order 7 2 ^ 3 , so much larger. Therefore, the approximation of 
f(x — £, y — 77-r, t) — f{x, y, t) in (|7^) by a second order Taylor expansion should still be correct for values of |£| and 
I r\ I up to C7 2 / 3 , with c some small positive constant. In addition it can be estimated that the contributions to the 
Boltzmann equation from |£|- and |?y|-values outside the cut-off are smaller by at least a factor of order (In i) 1 / 3 than 
the terms kept in the Fokker-Planck approximation. In this way, it is justified to introduce a cut-off of the form fls^), 
with 8 — |, which establishes the self-consistency of the scheme. 

The numerical computation shows that, indeed, the maximal Lyapunov exponent is of the form (|8^) for 7 <C 1, as 
seen in Fig. ||. The logarithmic correction turns out to be very small. We further observed that the second largest 
Lyapunov exponent also scales as (|39"|). 

Remark: The logarithmic divergences of the second moments of the distribution V are an anomaly of the hard-ball 
potential. For smooth potentials these moments are well-defined. 



B. The full fluid dynamics 

We now consider a fluid of small hard disks of radius a = h and mass m = 1 with a tracer disk of radius A and 
mass M. The temperature is always T = 1. In order to find the conditions under which the tracer particle dominates 
the Lyapunov spectrum we proceed as follows. 

We notice that the limit a — * is equivalent to the conditions A ^> a and na d <C 1. Therefore, we consider a 
sequence of systems of lower and lower density n. In the limit n — > 0, three types of tracer disks are considered: 



1G 



1 

2 " 


M = 


1, 


5, 


M = 


100, 


1 

2v^' 


M = 


10. 



(90) 

(ii) A = 5, M= 100, (91) 

(iii) A = —= . M = 10. (92) 

The specific purpose of the last choice is to consider sequences of Brownian particles of increasing sizes, but all made 
of the same material, so the ratio of M to nm(a + A) 2 remains constant. In Fig. the maximal Lyapunov exponent 
is plotted versus the density. 

Case (i) is the reference situation in which all disks are identical and we recover the pure fluid behavior ([32|), as 
expected. Indeed, a best fit of the numerical data to a linear combination of n Inn and n yields 

A f = -10.9 n Inn- 19.0 n, (93) 



with the coefficient of the leading term in agreement with the value predicted by Eq. (p2[): — 2y / 7rtj(41) = —10.89. 

Both cases (ii) and (iii) are examples of Brownian motion, as both A/ a and M/m are 3> 1. In case (ii) the parameter 
7 changes from ^> 1 for the larger of the density values considered to -C 1 for the lowest densities. In case (iii) 7 
remains almost fixed at a value of 4/11. The values of A/a and M/m chosen in case (ii) are such that the fluid always 
dominates the maximal Lyapunov exponent. This is clearly confirmed by the simulation results, which show that 
Ai ~ Af. 

In case (iii), in contrast, the tracer particle dominates the Lyapunov spectrum at low densities as observed in Fig. 
[| Indeed, because A — l/(2y/n), we infer from both Eq. fl57j ) and fl89|) that 

A t ~ u t ~ Vn (94) 

while the fluid Lyapunov exponent takes the values following from ( |93| ) so that Af <C A t for n — > 0. This explains the 
behavior observed in Fig. |9 We notice that case (iii) is intermediate between the two regimes, 7 > 1 respectively 



7 -C 1 studied in subsection [V A , because the mean free path of the tracer is only slightly smaller than the tracer 
radius. Indeed, in case (iii) we have £ t — 0.3 rT^ < A + a ~ 0.5 n - '. Therefore neither ( |57|) nor (|8|) strictly applies, 
but the scaling of the maximal Lyapunov exponent as y/n remains valid. 

However, with a system of one tracer disk among about 40 small disks, we are still in a situation very different from 
the Raylcigh flight because of the periodic boundary conditions and the boundedness of the domain. Even if the mean 
free path of the fluid particles among themselves is much larger than the size of the system Vf/vff 3> L Xl L y the fluid 
particles have recollisions on the tracer particle because of the periodic boundary conditions. When (A/2a) d ~ 1 Nf 
the dynamics of the fluid particles is dominated by recollisions with the tracer particle. We are then in a situation, 
similar to the Sinai billiard in which a point particle (such as a fluid particle) collides on a large disk (i.e., the tracer 
particle) in a system with periodic boundaries. This leads to two effects. 

The first effect is that, in the limit where the tracer particle is very massive, the fluid particles become almost 
independent of each other. Therefore the Lyapunov spectrum becomes similar to Nf copies of the Lyapunov spectrum 
(+As, 0, 0, —As) of the Sinai billiard. In the limit M — > 00, we should thus expect a Lyapunov spectrum of the form 
(+Ag, +As, 0, 0, —As, —As) in d = 2. This tendency is indeed observed in Fig. [io| which compares two systems 

JV f 2Ni N f 

of different sizes with an identical tracer particle and at the same density. We observe in Fig. |Io| that indeed the 
positive Lyapunov exponents roughly separate into two equally populated families. This tendency is stronger for the 
40-particle system than for the 80-particle one, as was to be expected in view of the preceding arguments. In the 
limit L x ,L y — > 00, we notice that the Lyapunov exponent of the Sinai billiard vanishes as As ~ (vf/L x>y )la.(L Xt y/A) 
so that this effect tends to disappear. 
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The second effect due to the periodic boundary conditions is more important for our present considerations than 
the previous one. If the recollisions of a fluid particle with the tracer particle become more frequent than collisions 
with other fluid particles these may have an important, or even dominant effect on the growth of the perturbations on 
the tracer coordinates. As a consequence, the maximal Lyapunov exponent of the system may be different from the 
Ray leigh- flight value even though the tracer particle dominates and controls it. This effect may even remove the gap 
in the Lyapunov spectrum, as seen in Fig. |ld. Notice that a small but finite density of large tracer particles will have 
exactly the same type of effects on the Lyapunov spectrum of a large or even infinite system as periodic boundary 
conditions in the case of a single tracer particle. 

In the limit where the system becomes infinite, we expect that both effects disappear and that a gap does appear in 
the Lyapunov spectrum. This is indeed the case, as observed in Fig. ^ showing the five largest Lyapunov exponents 
for systems with the same density n — 10~ 8 and temperature T = 1, the same tracer particle of radius A = 5000 and 
mass M = 10, but an increasing total size and so an increasing number N{ of fluid particles. Under these conditions, 
the limiting Ray leigh flight would have the following two positive Lyapunov exponents, obtained by Monte-Carlo 
simulation: 

d = 2, n= 10~ 8 , T = 1 , A = 5000 , M = 10 : 

A t = (3.3 ±0.1) x 10~ 5 , A{ = (0.6 ±0.1) x 10~ 5 (95) 

For the full dynamics simulated by molecular dynamics, we observe in Fig. |ll| that the maximal Lyapunov exponent 
decreases as N{ — > oo to a value close to the maximal Lyapunov exponent At of the Rayleigh flight. The next exponents 
decrease faster, creating a gap in the Lyapunov spectrum. We notice that, for the systems we studied, the second 
exponent has not yet converged to the Rayleigh-flight value X' t and is therefore not separated from the third and next 
exponents, but we expect this to occur for Nf large enough. 

In the Rayleigh-flight limit, a — > 0, combined with the large system limit N — > oo, we therefore find a Lyapunov 
spectrum that is dominated by the dynamical instability of the tracer particle: 

Ai ~ At > A 2 ~ A; > A 3 > A 4 > (96) 

with the formation of a gap, as in the case of the Lorentz-gas limit, but here because of a very different mechanism. 



V. CONCLUSIONS 



In this paper, we have studied a system of hard balls in elastic collisions (disks in d = 2 and spheres in d — 3) and we 
have shown that the tracer particle dominates the Lyapunov spectrum in the Lorentz-gas and in the Rayleigh-flight 
limit. 

In the Lorentz-gas limit, the tracer particle is lighter and moves faster than the fluid particles. The tracer particle 
therefore has a higher collision frequency than the other particles. Since the maximal Lyapunov exponent is propor- 
tional to the collision frequency, a gap appears in the Lyapunov spectrum between the largest Lyapunov exponents, 
which are associated with the tracer particle, and the rest of the spectrum. In d = 2, there is one such positive 
Lyapunov exponent associated with the tracer and, in d — 3, there are two such exponents. These largest Lyapunov 
exponents take values very close to the Lorentz-gas values previously obtained by Van Beijeren, Dorfman, and Latz 
H,^), as confirmed by direct numerical computation. 

The other limit in which the tracer particle dominates the Lyapunov spectrum is the Rayleigh-flight limit. In this 
limit, the radius of the fluid particles tends to zero or, equivalently, the density of the fluid particles vanishes while 
the radius of the tracer is much larger than the radius of the fluid particles. In an infinite system, the only collisions 
would occur between the fluid particles and the lone tracer particle. We have shown the remarkable result that, even 
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in this limit where the fluid is ideal and composed of non-interacting particles, the tracer particle may have positive 
Lyapunov exponents. We obtained formulas for the dependence of the maximal Lyapunov exponent on the parameters 
of the system in two different regimes. In the first regime, where the mean free path of the tracer particle is so much 
larger than its radius that perturbations of its coordinates are multiplied by large factors at almost all collisions, 
we find a behavior again very similar to that of the Lorentz gas. In the other regime, where the effect of a single 
collision on average is very small, we found that the maximal Lyapunov exponent scales as the product of the tracer 
particle collision frequency and the two-third power of the mass ratio m/(m + M), up to logarithmic corrections. 
These logarithmic corrections are special for the hard-ball potential and do not occur for more realistic interactions. 
Moreover, we have shown that the positive Lyapunov exponents determined by the tracer dynamics may dominate the 
Lyapunov spectrum of the fully interacting system under conditions approaching the Ray leigh- flight limit, provided 
the density of tracer particles, or alternatively the ratio between system size and tracer radius, remains finite. Again, 
a gap appears in the Lyapunov spectrum between the largest Lyapunov exponents associated with the tracer particle 
and the rest of the spectrum. 

The different tracer-dominated regimes are depicted in Fig. [jj] as a function of the mass ratio M/m and the 
radius ratio A/a for two different densities of a two-dimensional fluid. The diagrams are qualitatively similar for a 
three-dimensional fluid. Figure n3shows that the system is in a fluid-dominated regime if the tracer is very massive. 
Nevertheless, we observe in Fig. ]12|b that the tracer-dominated regimes extend toward larger masses at lower densities. 

A comment is here in order about Brownian motion. In typical Brownian motion conditions, the tracer particle is 
much more massive than the fluid particles and, moreover, the mean free path of the fluid particles among themselves 
is much shorter than the radius of the Brownian particle. Under such conditions the maximal Lyapunov exponent of 
the full system usually is essentially the same as the fluid Lyapunov exponent, as seen in Fig. [l^. Then the Brownian 
particle does not contribute significantly to the dynamical instability of the system and is a probe for the dynamics 
of the surrounding fluid. This is the case in typical Brownian-motion experiments. In order to observe the new 
effect of dominance of the dynamical instability by the Brownian particle one has to use a sufficiently rarefied gas as 
surrounding fluid, in order to approach to the Rayleigh-flight limit, as seen in Fig. [T^b. 
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FIG. 1. Lyapunov spectra of a 2d fluid of density n = 1CP 3 and temperature T = 1 composed of iVf = 40 hard disks of 
radius a = 1/2 and mass m = 1 with 7V t = 1 hard disk of radius A = 1/6 and mass: (a) M = 10; (b) M = 10~ 2 . The squares 
depict the positive exponents and the crosses are minus the negative exponents. 
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FIG. 2. Collision frequency u t of the tracer particle in a 2d fluid of density n — 10 and temperature T = 1 composed of 
N{ = 40 hard disks of radius a — 1/2 and mass m — 1 with JVt = 1 hard disk of radius A — 1/6 and varying mass M. The 
filled circles are the numerical data. The solid line is the prediction of Eq. ( ^o[ ) (see Table I). 
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FIG. 3. The five largest Lyapunov exponents of a 2d fluid of density n = 10 and temperature T = 1 composed of JVf = 40 
hard disks of radius a — 1/2 and mass m = 1 with JVt = 1 hard disk of radius A = 1/6 and varying mass M. The filled circles 
depict the largest Lyapunov exponent and the open circles the four next ones. The solid line is the prediction of Eq. (MW. 
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FIG. 4. The five largest Lyapunov exponents of a 2d fluid of density n = 10~ and temperature T = 1 composed of iVf hard 
disks of radius a = 1/2 and mass m — 1 with jV t = 1 hard disk of varying radius A and mass M = 10 -2 (JVf = 40 for 4 < 10, 
ATf = 39, 37, and 31 for the last three values of A > 10). The filled circles depict the largest Lyapunov exponent and the open 
circles the four next ones. The solid line is the prediction of Eq. (]4q). 
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FIG. 5. Lyapunov spectrum of a 3d fluid of density n = 10 -3 and temperature T = 1 composed of Nt = 48 hard spheres 
of radius a — 1/2 and mass m = 1 with N t — 1 hard sphere of radius .4 = 1/6 and mass M = 10 -2 . The squares depict the 
positive exponents and the crosses are minus the negative exponents. 
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FIG. 6. The five largest Lyapunov exponents of a 3d fluid of density n = 10 -3 and temperature T = 1 composed of Nt = 48 
hard spheres of radius a — 1/2 and mass m = 1 with JVt = 1 hard sphere of radius A — 1/6 and varying mass M. The filled 
circles depict the two largest Lyapunov exponents and the open circles the three next ones. The solid lines are the predictions 
of Eqs. @ and @. 
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FIG. 7. Typical geometry of two successive collisions of the tracer particle when: (a) i\ A + a; (b) l t <§C A + a. The 
trajectories are depicted in the frame where the fluid particle is at rest. 
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FIG. 8. Ratio of the maximal Lyapunov exponent At to the collision frequency v% versus the mass M for a tracer disk of 
radius A = 499.5 in Rayleigh flight in the regime y <C 1. The Lyapunov exponent is calculated by Monte-Carlo simulation by 
assuming that the tracer particle undergoes independent successive random collisions from fluid disks of radius a — 0.5 and 
mass m = 1 at density n = 10™ 8 and temperature T = 1. The straight line has the theoretically expected slope 2/3. 
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FIG. 9. maximal Lyapunov exponent of a 2d fluid of temperature T = 1 and varying density n containing Nf ~ 40 hard 
disks of radius a — 1/2 and mass m = 1 with one tracer disk of radius A and mass M in the three cases (^). In the reference 
case A = 1 and M = 1, the maximal Lyapunov exponent Ai is depicted by the crosses and the fit by the solid line. In the case 
A — 5 and M — 100, Ai is depicted by the open squares and the fit by the dashed line. In the case A — l/(2^n) and M — 10, 
Ai is depicted by the filled squares. The long-short dashed line is the fit Ai = 0.75-^71. 
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FIG. 10. Lyapunov spectra of a 2d fluid of temperature T = 1 and density n = 10 containing respectively Nf = 39 (large 
squares) and Nf = 83 (small squares) hard disks of radius a = 1/2 and mass m = 1 with one tracer disk of radius A — 5000 
and mass M — 10 as a function of the relative index i/(2N) of the Lyapunov exponents Aj. 
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FIG. 11. The five largest Lyapunov exponents of a 2d fluid of temperature T = 1 and density n = 10 containing a varying 
number Nf — 39, 83, 143, 239, 359, 503, 671, 863 of hard disks of radius a = 1/2 and mass m = 1 with one tracer disk of radius 
A — 5000 and mass M = 10 versus 1/iVf. The Rayleigh-flight value of the maximal Lyapunov exponent is A t ~ 3.3 10~ 5 . 
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FIG. 12. Diagrams of the different regimes of the d = 2 systems composed of a tracer particle of radius A and mass M in a 
fluid of particles of radius a and mass m at density (a) n = and (b) n = 10~ 10 and temperature T = 1. The area above 

the solid line shows the fluid-dominated regime with Af > At. The areas below the solid line show the tracer-dominated regimes 
with At > Af. The area on the lower left-hand side is the tracer-dominated regime with 7 > 1, while the area on the lower 
right-hand side is the tracer-dominated regime with 7 < 1. The Lorentz-gas limit is the part of the lower area where M <§C m. 
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